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1 Introduction 

MEG and EEG have the time resolution necessary 
to look at the visual system’s real-time response to 
alternating visual stimuli. Currently, investigators 
in the field of EEG have carried out the most 
complete study of the human brain's response to 
alternating visual stimuli. Some of the early work 
using steady-state visually evoked potentials 
(SSVEP) stressed the importance of controlling all 
physical properties of the stimulus, van der Tweel 
et al. [1] used a sinusoidally modulated white light 
source at many different modulation frequencies. 
They found that the sinewave shape was 
reproduced very well in the EEG at 9-15 Hz for 
diffuse large field illumination. Regardless of the 
purely sinusoidal input, 2nd harmonics were found 
in response to diffuse illumination between 5-9 Hz. 
Only one study measuring SSVEPs attempted 
rudimentary localization using an available EEG 
localization tool, VARETA [2]. This selective 
attention study employing SSVEPs (20-28 Hz 
range) found that contralateral occipital-parietal 
cortex was primarily responsible for the attention 
effects noted. Due to the current limitations in the 
localization techniques in EEG, these studies have 
not identified the generators of the oscillatory 
response itself. 

Okada et al. [3] conducted a comprehensive study 
comparing MEG and EEG responses. These results 
were consistent with the EEG findings in that the 
visual system responded well to stimuli at 3-30 Hz. 
Some preliminary localization of the alternating 
response was also performed assuming a single 
equivalent current dipole that localized the source 
to the occipital cortex. More recent time-varying 
visual MEG studies have either focused on moving 
stimuli [4,5,6,7] or on alternating stimuli at a rate of 
1 Hz (i.e. transitory responses) [8,9]. 

To our knowledge, this is the first study to employ 
automated multi-dipole spatio-temporal modeling 
techniques to determine the origins and the 
corresponding timecourses of the responses to 


alternating visual stimuli. This allowed for 
individual identification of sources responsible for 
the oscillatory behavior found in the magnetic field 
responses. 

2 Methods 

2.1 Subjects and Stimuli 

Five normal subjects (28-48 years; 2 female) were 
studied. Informed consent was obtained from all 
subjects prior to study involvement. Sixteen 
different visual stimuli were presented along the 
horizontal meridian to the subjects. The visual 
stimuli were back-projected through a small 
window from an Epson LP735 projector onto a 
rear-projection lenticular screen placed in front of 
the subject. Four different frequencies (0, 3, 8.7 and 
14 Hz) of alternating black and white circular 
sinusoids with a luminance matched background 
were presented to left and right foveal (eccentricity 
2.3°, visual angle 1.8°, spatial frequency 2.6/cm) 
and peripheral locations (eccentricity 24°, visual 
angle 7°, spatial frequency 0.6/cm). 

2.2 Data Collection 

The MEG data was co-localized to the anatomical 
MRI of each subject using information obtained 
from a Polhemus head position device. The 
magnetic field responses were measured with a 
whole-head 122-channel Picker/Neuromag 
biomagnetometer. The conditions were randomized 
to minimize drift away from the central fixation 
point. The stimulus duration was one second, 
thereby presenting a sufficient number of cycles of 
all frequencies to appear periodic. The ISI was 1.7 
± 0.2 seconds. The responses were on-line signal 
averaged (N = 200) at a digitization rate of 300 Hz 
with a band-pass filter of 0.03-100 Hz. An interval 
including 100 ms pre-stimulus, baseline data and 
1500 ms post-stimulus data was collected. 



2.3 Analysis 

All of the data obtained were preprocessed using 
MEGAN (MEG ANalysis), a MEG data analysis 
tool developed by the biophysics group at Los 
Alamos National Laboratory. The data were filtered 
with a band-stop 60 Hz filter. All bad channels 
were removed (i.e. channels with flux jumps, or 
large amplitude responses found in only one 
channel). 

An initial analysis of an early interval of time (60- 
180 ms post-stimulus) was performed to verify the 
quality of the MEG data. Once the quality of the 
visual response was established, the waveforms 
were inspected for their frequency following (ff) 
characteristics using MEGAN. The ff portion of the 
response occurred in the 250-1000ms post-stimulus 
time interval. The waveforms were rated ff if any 
channels contained significant power at the 
harmonics of the driving frequency. Any power 
value greater than the average 0 Hz power plus 
three standard deviations was considered 
significant. 

The data were analyzed using a multi-dipole spatio- 
temporal modeling technique [10] that performs a 
downhill simplex search for hundreds to thousands 
of fits to the data using randomly selected starting 
parameters. To determine the number of sources to 
be modeled, singular value decomposition (SVD) 
was performed on each data set for the time interval 
of interest. The dipole fitting procedure was 
performed for several different numbers of dipoles, 
or model orders, around the estimated model order 
obtained with the SVD. The most appropriate 
model order was found by comparing the goodness- 
of-fit and comparing the reduced-chi-square values 
of the different model orders [11,12]. 

Once the best solutions were obtained for each 
condition and for each of the five subjects, the 
timecourses of each dipole location were rated for 
their ff characteristics. The power at the first and 
second harmonics was found, to determine the 
extent of ff involvement using a normalized power 
spectrum function in Matlab. Any timecourse with 
power greater than the 0 Hz mean power plus 3 
sigma at the fl or f2 frequencies was labeled ff. The 
ff results were then collated across conditions and 
subjects. The significance values found in Table 1 
were determined using chi-square statistics for 
comparing frequency of occurrence. 

3 Results 

Upon initial inspection of the raw waveforms, some 
conditions showed obvious oscillatory behavior as 
seen in Figure 1. The frequencies apparent in the 


2 nd , 3 rd , and 4 th panels are approximately 6, 17.4 and 
14 Hz, respectively. These frequencies correspond 
to the 2 nd harmonic (f2) of the 3 and 8.7 Hz 
stimulus and the 1 st harmonic (fl) of the 14 Hz 
stimulus. Although oscillatory behavior was not 
obvious in many of the conditions, inspection of 
individual channels and the corresponding power 
spectra confirmed ff in most of the conditions 
across subjects. There was no clear increase in 
power above the second harmonic for any of the 
frequencies. The fl and f2 pattern seen in the 
waveforms was consistent across subjects. 
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Figure 1: Examples of raw waveforms from the four 
frequency conditions. The responses from all 122 
channels in each plot are overlaid. No clear 
oscillatory response is seen in the 0Hz condition; 
whereas, 6 and 17.4 Hz oscillations are seen in the 
3 and 6 Hz conditions, respectively. The prominent 
oscillations in the 14 Hz condition are oscillating at 
MHz. 

Figure 2 shows the results of a multi-dipole fit to 
the late interval of time in one condition from one 
subject. The appropriateness of the 6-dipole model 
can be seen by noting the similarity between the 
raw waveforms (Empirical) and the modeled 
waveforms (forward). In addition, one can see 
spatially distinct areas in the three-dimensional plot 
as well as temporally distinct timecourses. Notice 
that some of the timecourses exhibit oscillatory 
behavior and some do not. 

Examples of ff timecourses can be seen in Figure 3. 
Corresponding examples of the areas exhibiting ff 
behavior are shown in Figure 4. Again, the fl and 
f2 results are consistent with the results of the 
waveforms, although in addition there were areas 
that responded at 28 Hz, the f2 of the 14Hz 
stimulus. 
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Figure 2: a) The four panels show the results of a 6 
dipole fit for the 250-800 ms time interval for the 3 
Hz left central field stimulation in subject MH. The 
residual is the difference between the Empirical 
and forward waveforms, b) Shows a 3D plot of the 
dipole locations of the 5 best fits of the multi-start 
results. The tight clusters of asterisks from the five 
different fits suggest a global minimum was 
obtained, c) Shows the timecourses of the six dipole 
locations. Red - L. Inferior Temporal, Green - R. 
STS, Blue - R. Lateral Occipital, Magenta - R. MT, 
Aqua - L. Frontal, Yellow - L. Posterior Cingulate. 
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Figure 3: Examples of ff timecourses from different 
subjects. The top two rows show ff timecourses 
oscillating at 6 and 17.4 Hz in response to the 3 
and 8.7 Hz conditions, respectively. The bottom two 
rows show timecourses oscillating at 14 and 28 Hz 


in response to the 14 Hz condition. They-axis is the 
strength of the dipole in units of nA-m. Six of these 
timecourses were associated with V1/V2 activity. 
Two of the timecourses were associated with 
activity from area MT (2 lower right). 
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Figure 4; Locations found to have a ff response. 
The different areas are represented by different 
subjects. MT - medial temporal-occipital junction, 
POS - parietal-occipital sulcus, TPJ - temporal- 
parietal junction. 


Table 1: Prevalence of ff areas compared to the 
prevalence of all areas active in the late time 
interval (250-800 ms). ___ 



V1/V2 

Parietal 

MT 

POS 

ff (N=84) 

36%* 

7% 

13%** 

8% 

all (N=345) 

22% * 

12% 

7%** 

13% 


FEF 

Cing. 

TPJ 

*P=0.05 

**P=0.1 

ff (N=84) 

6%* 

15% 

10% 

all (N=345) 

12%* 

13% 

12% 

















































































Table 1 shows the prevalence of the ff response in 
the different areas active during the late interval of 
time. This is compared to the prevalence of all 
areas active in the late time interval. Clearly, the 
V1/V2 area is primarily responsible for the ff 
response. Six of the seven areas were represented 
by 4/5 subjects in the ff timecourses. The temporal 
parietal junction area was only represented by 3/5 
subjects. Putative area MT and cingulate cortex is 
also quite well represented in the ff timecourses. 

4 Discussion 

The variability found in the raw waveforms across 
subjects is consistent with the literature [13,14]. 
Although there is considerable variability in the 
literature concerning the visual system’s response 
to alternating stimuli at the various harmonics, 
regardless of experimental design, the appearance 
of the first harmonic at high frequencies is 
consistent with the overall trend. 

The main result of this study was that the frequency 
following response was primarily generated by 
primary visual cortex as suggested by [3]. 
However, the presence of ff timecourses in higher 
visual areas suggests that these areas are also 
involved. It is likely that the ff nature is identified 
less often in the higher order areas due to a weaker 
response - leading to less reliable modeling. In 
addition, middle temporal occipital cortex (MT) is 
sensitive to drifting stimuli but somewhat less 
sensitive to alternating stimuli in monkeys. Similar 
results have been found in PET and fMRI studies, 
as well as a few MEG studies [4,5], suggesting that 
the type of stimulus may also lead to weak ff in 
area MT. However, it is interesting to note that half 
of the MT sources found in the late activity were 
labeled ff compared to 40% for V1/V2. Cingulate 
cortex has most often been implicated in 
maintaining attention. With alternating stimuli, it is 
likely that cingulate cortex is active to monitor the 
incoming stimuli. These results suggest that the raw 
frequency information is relayed to the higher 
visual areas. 
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